Matter-wave analog of an optical random laser 
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The accumulation of atoms in the lowest energy level of a trap and the subsequent out-coupling of 
these atoms is a realization of a matter- wave analog of a conventional optical laser. Optical random 
lasers require materials that provide optical gain but, contrary to conventional lasers, the modes 
are determined by multiple scattering and not a cavity. We show that a Bose-Einstein condensate 
can be loaded in a spatially correlated disorder potential prepared in such a way that the Anderson 
localization phenomenon operates as a band-pass filter. A multiple scattering process selects atoms 
with certain momenta and determines laser modes which represents a matter-wave analog of an 
optical random laser. 

PACS numbers: 03.75.Pp, 42.55.Zz 



Conventional optical lasers require two ingredients: 
material that provides optical gain and an optical cavity 
responsible for coherent feedback and selection of reso- 
nant laser modes. However, it is also possible to achieve 
laser action without the optical cavity provided the gain 
material is an active medium with disorder [ij. Forty 
years ago Letokhov analyzed a light diffusion process 
with amplification and predicted that gain could over- 
come loss if the volume of a system exceeded a criti- 
cal value 0. Random lasing (i.e. light amplification in 
disordered gain media), achieved in a laboratory in the 
1990's, attracts much experimental attention and offers 
possibilities for interesting applications Theoret- 
ical understanding of this phenomenon is still imperfect. 
Although the Letokhov model of diffusion with gain is 
useful in predicting certain properties of random lasers, 
it neglects coherent phenomena. There are various theo- 
retical models of random lasing but it is widely accepted 
that interference in a multiple scattering process deter- 
mines the spatial and spectral mode structure of a ran- 
dom laser [1|. 

Bose-Einstein condensation (BEC) of dilute atomic 
gases is a macroscopic accumulation of atoms in the low- 
est energy level of a trap when the temperature of the 
gas decreases Q. This tendency of occupying a single 
state through the mechanism of stimulated scattering of 
bosons is an analog of mode selection in optical lasers due 
to the stimulated emission of photons. Gradual release 
of atoms from a trapped BEC allows for the realization 
of a matter-wave analog of a conventional optical laser 
[§-14 1 . The atom trap is an analog of the optical cav- 
ity. The lowest mode of the trap is a counterpart of an 
optical resonant mode. In conventional optical lasers the 
output coupler is usually a partially transmitting mirror. 
In atom lasers it involves, for example, a change of the 
internal state of the atom by means of a radio-frequency 
transition. 

In the present letter we propose the realization of a 
matter-wave analog of an optical random laser. Suppose 



the BEC of a dilute atomic gas has been achieved in a 
trapping potential. That is, we begin with the accumu- 
lation of atoms in a single mode of the resonator (i.e. the 
lowest cigcnstate of the trap). Then, let us turn off the 
trap and turn on a weak disorder potential. Starting with 
a BEC we have a guarantee that the disorder medium is 
pumped with coherent matter- waves. We would like to 
raise the question of whether it is possible to prepare 
spatially correlated disorder potential in such a way that 
narrow peaks can be observed in the spectrum of atoms 
that are able to leave the area of the disorder potential? 
In other words: if the multiple scattering of atoms in a 
disorder medium can lead to a selective spectral emission 
of matter- waves from the medium? 

In cold atom physics a disorder potential can be re- 
alized by means of an optical speckle potential [l^ [3 ■ 
Transmission of coherent light through a diffusing plate 
leads to a random intensity pattern in the far field. 
Atoms experience the presence of the radiation as an ex- 
ternal potential V^(r) oc x|-E'(r)p proportional to the in- 
tensity of the light field E{r) and atomic polarizabitlity x 
whose sign depends on the detuning of the light frequency 
from the atomic resonance. Diffraction from the diffusive 
plate onto the location of atoms determines correlation 
functions of the speckle potential. We assume that the 
origin of the energy is shifted so that ^(r) = where 
the overbar denotes an ensemble average over disorder 
realizations. Standard deviation of the speckle potential 
Vo measures the strength of the disorder. 

Let us begin with a one-dimensional (ID) problem. In 
a weak disorder potential atoms with /c-niomentum un- 
dergo multiple scattering, diffusive motion and finally lo- 
calize with an exponentially decaying density profile due 
to the Anderson localization process, provided that the 
system size exceeds the localization length [l7l - [l9j . Tak- 
ing the Born approximation to the second order in the 
potential strength, the inverse of the localization length 
is = [mVo/h^kf'P{2k), where the Fourier trans- 

form of the pair correlation function of the speckle po- 
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and 7(fc) = J dz^{z)e is the Fourier trans- 
form of the complex degree of coherence 7(2:) — 
E*{z + z')E{z')/\E{z)\^ = / dyA{v)e^^yl°'/ J dyA{y). 
A{y) describes the aperture of the optics and a is a 
constant dependant on the wavelength of the laser ra- 
diation and the distance of the diffusive plate from the 
atomic trap. In Ref. 2lJ, where the experimental real- 
ization of the Anderson localization of matter-waves is 
reported (sec also [l^]), a simple Heavisidc step function 
A{z) = Q{R — \z\) describes the aperture. The corre- 
sponding 7(fc) = 7r(TflO(l — |A:(Tii|), where cte. = a/R 
is the correlation length of the speckle potential. Con- 
sequently the power spectrum ([T]) decreases linearly and 
becomes zero for > 2/aii- Thus, the Born approxima- 
tion predicts an effective mobility edge at |fc| = i.e. 
atoms with larger momenta do not localize [21I [2^ (ac- 
tually higher order calculations show they do local- 
ize but with very large localization lengths, much larger 
than the system size in the experiment). Hence, neglect- 
ing atom interactions, if the width of the initial atom 
momentum distribution exceeds the mobility edge par- 
ticles at the tail of the distribution avoid the Anderson 
localization and may leave the disorder area. 



Let us modify the experiment reported in Ref. [21| by 
introducing an obstacle at the center of the diffusive plate 
so that the aperture is now described by A{z) = 0(i? — 
\z\) — Q{p — \z\) where p < R. It implies that 
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where o-p = a/ p. If the size of the obstacle p > i?/3, in- 
terference of light passing through such a double-slit dif- 
fusive plate creates a peculiar speckle potential. That is, 
the power spectrum ([1]) disappears for |fc| > 2/ an as pre- 
viously but it is also zero for — — < \k\ < —. Thus, 

according to the Born approximation there is a momen- 
tum interval where the localization length diverges. It 
implies that the Anderson localization process is able to 
operate as a band-pass filter letting particles with spe- 
cific momenta leave the region of the disorder. Detection 
of escaping atoms should reveal a peak in the momen- 
tum spectrum corresponding to the interval where the 
localization length diverges. 

Introducing two (or more) obstacles in the diffusive 
plate we can increase the number of momentum intervals 
with diverging lioc- In Fig. [T] we present examples of the 
speckle potentials in the single obstacle case and the case 
of two obstacles located symmetrically around the plate 
center. In the latter case the aperture is described by 
Aiz) = e{R- \z\) - e(p- 1^1) + e(C- |z|) where C < P 
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Figure 1: (Color online) Examples of the speckle potential 
(top panels) and the corresponding localization length (bot- 
tom panels) obtained within the Born approximation (dashed 
black curves) and numerically in the transfer-matrix calcula- 
tions (solid red curves). Panels (a) and (c) show the results 
for the single obstacle in the diffusive plate where an = 0.066 
(0.31 ^m), Ga/Gp = 0.4 and Vb = 3.5. Panels (b) and (d) cor- 
respond to the case of two obstacles with the same gr and Vo 
but gr/gp — 0.7 and gr/g(^ = 0.1. The results are shown for 
rubidium-87 atoms. Red detuning of the laser radiation from 
the atomic resonance is assumed. All values are presented in 
the harmonic oscillator units, i.e. energy Eq — hoj and length 
lo ~ ^/h/mw where oj/2tv = 5.4 Hz. 
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with — a/C,. The figure also presents localization 
lengths obtained numerically in the transfer-matrix cal- 
culation [iJl that confirms the Born predictions. 

To simulate an experiment, we follow the parameters 
used in Ref. 2l| where Anderson localization of matter- 
waves has been observed. We assume that a BEC of 
N ~ 1.7 • 10* rubidium-87 atoms is initially prepared 
in a quasi- ID harmonic trap with longitudal and trans- 
verse frequencies 10/2-11 = 5.4 Hz and a;_L/27r = 70 Hz, 
respectively. In the following we adopt the harmonic os- 
cillator units: Eq = Huj, Iq = y/h/riiu and to ~ 1/w for 
energy, length and time, respectively. When the trap- 
ping potential is turned off and the speckle potential is 
turned on the expansion of the atomic cloud is initially 
dominated by the particle interactions until the density 
drops significantly and the atoms start feeling only the 
disorder potential. This initial stage of the gas expan- 
sion sets the momentum distribution of the atoms which 
may be approximated by an inverted parabola with an 
upper cut-off k,nax = 2y/JI = 12.7 where p, is the initial 



chemical potential of the system [23|, ]2 



The disorder potentials we choose are attainable in the 
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Figure 2: (Color online) Momentum distributions of the 
atoms localized in the disorder potential (dashed black lines) 
and the atoms which escaped from the disorder area (solid 
red lines). Panels (a) and (c) show the results for a single ex- 
perimental realization of the disorder with parameters as in 
Fig. [TJi,c while panels (b) and (d) are related to the param- 
eters as in Fig. [T]D,d. Panel (a) corresponds to the evolution 
time t = 100 (2.9 s), panel (b) to t = 70 (2 s), panels (c) and 
(d) tot = 200 (5.7 s). Fraction of atoms that escaped the dis- 
order region is about 9% in (a) and (b), and 20% in (c) and 
(d). In order to take into account experimental resolution all 
data have been convoluted with Gaussian of Afc = 0.3 width. 
All values are presented in the harmonic oscillator units, i.e. 
energy Eq — huj, length lo — ^h/mu! and time to = l/co 
where cj/27r — 5.4 Hz. 



experiment reported in Ref. (2l|, i.e. they extend 862 
units (4 mm) along the z direction with the correlation 
length cTfl = 0.066 (0.31 /im). We consider the potentials 
obtained by introducing one or two obstacles in the dif- 
fusive plate that are presented in Fig. [TJ If the atomic 
cloud starts at the center of the disorder potentials and 
if the cut-off of the momentum distribution is kmax ^ 13 
then we may expect that time evolution leads to emis- 
sion of atoms with |fc| w 5.5 in the single obstacle case 
(cf. Fig. [T^,c) and with |fc| ss 9 in the case of two obsta- 
cles (cf. Fig.[Tl3,d). In the latter case we may expect also 
a small leakage of atoms with |fc| w 3.5. For |fc| w 3.5 
the localization length shown in Fig[T}i reaches locally a 
maximum value of hoc ~ 120. Before the particle interac- 
tions become negligible the gas spreads over a significant 
range of the disorder region. Therefore the Anderson lo- 
calization is not able to diminish completely the leakage 
of atoms with |fc| sa 3.5. 

Starting with the ground state of the stationary Gross- 
Pitaevskii equation [9l4l5| in the presence of the harmonic 
trap, we integrate the time-dependent Gross-Pitaevskii 
equation when the trap is turned off and a disorder is 
turned on. Figure [2] shows momentum distributions of 
atoms that escaped from the disorder region and those 
that are localized for the disorder potentials correspond- 
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Figure 3: (Color online) The Boltzmann transport mean- 
free path (a) and the localization length (b) for atoms in 
the 2D speckle potential created by transmission of a laser 
beam through the circularly shaped diffusive plate with the 
obstacle in the form of a ring, i.e. the aperture of the op- 
tics is described by Q{R - |r|) - Q{p - |r|) -I- e(C - |r|) with 
ffl/fp = 0.99 and an/ac^ = 0.15. The potential strength cor- 
responds to »7 = 0.15. The quantities presented in the figure 
are dimensionless. 



ing to Fig. [T] at different moments in time. The expected 
selective spectral emissions of atoms are apparent in the 
figure. Interestingly in Fig. [2}i, i.e. for longer evolution 
time, small peaks around |fc| w 3.5 become visible. 

Finally let us consider a possibility of the realization 
of an atom analog of an optical random laser in 2D. The 
Boltzmann transport mean- free path Ib is the charac- 
teristic spatial scale beyond which memory of the ini- 
tial direction of the particle momentum is lost. In 2D 
hoc = ^se'^'^'^/^ and thus the localization length is much 
larger than Is which is in contrast to the ID case where 
these two quantities are nearly identical koc ~ '21b 
For the circularly shaped diffusing plate with the radius 
R the classical transport mean- free path [2^ [2^ , to the 
second order in the potential strength, reads 



(4) 



where rj = Vq/Ect is the ratio of the potential strength 
and correlation energy Eg. = h^/{maj^) with aji = a/R. 
The power spectrum 'P{k) of the optical speckle poten- 
tial disappears for k > 2 /an. Nevertheless, the Ib (and 
consequently also hoc) is always finite. In the bulk 2D sys- 
tem, an initially prepared atomic wave-packet follows a 
diffusive motion at short time but eventually the dynam- 
ics slow down and freeze due to the Anderson localization 
process [25l - [27t . 

By introducing obstacles in the diffusive plate we are 
able to shape the power spectrum of the speckle poten- 
tial. On one hand the fact that V{k) may disappear at 
certain momentum intervals does not mean divergence 
of the corresponding transport mean-free path On 



4 



the other hand any non-monotonic behaviour of Isik) is 
dramatically amplified in the behaviour of lioc(k) because 
the localization length is an exponential function oi Ib- 
In Fig. [3] we present an example related to the obstacle 
in the form of a ring, i.e. the aperture of the optics is 
described by ^(r) = e(i? - |r|) - e(p - |r|) + e(C - |r|). 
At kan « 0.4 both Is and lioc shows a maximum. How- 
ever, while the transport mean-free path changes by only 
a few in the neighboring region the localization length 
changes by four orders of magnitude. If the width of 
the momentum distribution of a BEC loaded in such a 
disorder potential is smaller than Q.Q/aji and the radius 
of the disorder medium is greater than IQ^an but less 
than lO^CTii; the multiple scattering process leads to an 
isotropic emission of atoms with k w QA/aji. 

We have outlined a proposal for the realization of a 
matter-wave analog of an optical random laser. Spatially 
correlated disorder potential for atoms with a peculiar 
pair correlation function can be created by transmitting 
a laser beam through a diffusive plate with obstacles. 
The resulting Anderson localization length reveals non- 
monotonic behaviour as a function of particle momen- 
tum. It allows for filtering momenta of particles that 
leave the area of the disorder, if the size of the disor- 
der medium is suitably chosen. The disorder medium is 
assumed to be initially loaded with a BEC which guaran- 
tees that the matter-waves emitted from the medium are 
coherent. We have restricted ourselves to the ID and 2D 
cases but the atom analog of an optical random laser can 
be also anticipated in 3D. In 3D the loffe-Regel criterion 
discriminates between waves that are Anderson localized 
{klB < 1) or not Thus, a spatially correlated dis- 



order potential for which the loffe-Regel criterion is not 
fulfilled for specific momenta should allow for selective 
emission of matter- waves in 3D. 

Our proposal is directly applicable to atomic matter- 
wave experiments. From the point of view of the optical 
random lasers our analysis is not complete because it 
is restricted to passive random materials without gain. 
There is an interesting question whether a disorder with 
properties similar to those analyzed here play a role in 
optical random lasers and which modes are important 
when the gain is included in a system. 

The non-monotonic behaviour of the localization 
length results in the appearance of a multiple effective 
mobility edge if a disorder system is finite. Wave trans- 
port is then unusual and interesting on its own. A shallow 
non-monotonical bahaviour of the Anderson localization 
length versus energy has been observed also in a classical 
wave system, see Ref. j29j . 
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Note added: After submission of this article, we be- 
came aware of a related theoretical study [s^ . 



[1] 
[2] 

[3] 
[4] 

[5] 

[6] 

[7] 
[8] 
[9] 
[10] 

[11] 
[12] 
[13] 
[14] 
[15] 
[16] 
[17] 
[18] 

[19] 



[20] 



[21] 
[22] 
[23] 

[24] 

[25] 
[26] 
[27] 

[28] 
[29] 
[30] 



D. S. Wiersma, Nature Physics 4, 359 (2008). 

V. S. Letokhov, Zh. Eksp. Teor. Fiz. 53, 1442 (1967); 

Sov. Phys. JETP 26, 835 (1968). 

N. M. Lawandy et al, Nature 368, 436 (1994). 

D. S. Wiersma and A. Lagendijk, Phys. Rev. E 54, 4256 

(1996). 

H. Cao et al, Appl. Phys. Lett. 73, 3656 (1998); H. Cao 
et al, Phys. Rev. Lett. 82, 2278 (1999). 

S.V. Frolov et al, Phys. Rev. B 57, 9141 (1999); T.V. 
Shahbazyan et al, Phys. Rev. B 61, 13 266 (2000). 

D. S. Wiersma, Nature 406, 135 (2000). 

M. H. Anderson et al, Science 269, 198 (1995). 
M.-O. Mewes et al, Phys. Rev. Lett. 78, 582 (1997). 

E. W. Hagley et al, Science 283, 1706 (1999). 

I. Bloch et al, Phys. Rev. Lett. 82, 3008 (1999). 

G. Cennini et al, Phys. Rev. Lett. 91, 240408 (2003) 
W. Guerin et al, Phys. Rev. Lett. 97, 200402 (2006) 

A. Couvert et al, Europhys. Lett. 83, 50001 (2008). 

J. T. Schulte et al, Phys. Rev. Lett. 95, 170411 (2005). 

D Clement et al. New J. Phys. 8, 166 (2006). 

PW. Anderson, Phys. Rev. 109, 1492 (1958). 

P. A. Lee and T.V. Ramakrishnan, Rev. Mod. Phys. 57, 

287 (1985). 

B. van Tiggelen, in Diffuse Waves in Complex Media, 
edited by J.-P. Fouque, NATO Advanced Study Insti- 
tutes, Ser. C, Vol. 531 (Kluwer, Dordrecht, 1999). 

I. M. Lifshits, S. Gredeskul, and L. A. Pastur, Introduc- 
tion to the Theory of Disordered Systems (Wiley, New 
York, 1988). 

J. Billy et al. Nature 453, 891 (2008). 

G. Roati et al. Nature 453, 895 (2008). 

L. Sanchez-Palencia et al, Phys. Rev. Lett. 98, 210401 

(2007) 

E. Gurevich and O. Kenneth, Phys. Rev. A 79, 063617 
(2009); P. Lugan et al, Phys. Rev. A 80, 023605 (2009). 
R. C. Kuhn et al, Phys. Rev. Lett. 95, 250403 (2005) 

C. Miniatura et al, Eur. Phys. J. B 68, 353 (2009). 

M. Robert-de-Saint- Vincent et al, Phys. Rev. Lett. 104, 
220602 (2010). 

R. G. Kuhn et al. New J. Phys. 9, 161 (2007). 

D. Capeta et al, Phys. Rev. A 84, 011801(R) (2011). 
M. Piraud, A. Aspect, L. Sanchez-Palencia, 
iarXiv:1104.2314. 



